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Abstract 

Aircraft  images  were  collected  near  Lee  Stocking  Island  (LSI),  Bahamas,  with  wavelike  features  for  bright  sand 
bottoms  during  times  when  solar  zenith  angles  were  large.  The  image  contrast  between  leading  and  trailing  wave 
facets  approached  a  10-15%  difference  because  of  algae  accumulations  in  wave  troughs  or  topographic  variations 
of  the  bottom.  Reflectance  contrast  for  blue  light  was  greater  than  for  red  and  green  wavelengths  when  algae  or 
detritus  was  present  in  the  ti'oughs.  However,  the  contrast  at  green  and  red  wavelengths  was  greater  than  at  blue 
wavelengths  when  caused  by  the  interplay  between  bottom  topography  and  oblique  illumination.  A  three-dimen¬ 
sional  backwards  Monte  Carlo  (BMC)  model  was  developed  to  evaluate  the  effect  of  oblique  illumination  on 
wavelike  topographic  features  for  various  values  of  water  clarity  and  bottom  albedo.  An  inverse  optical  modeling 
approach,  previously  developed  for  flat,  horizontally  homogeneous  bottoms,  was  applied  to  the  BMC  results.  Bathy¬ 
metric  estimates  for  bright  facets  tilted  10°  toward  the  sun  were  slightly  smaller  than  actual  depths,  whereas  shaded 
facet  depth  estimates  were  too  high  by  about  5%.  Larger  errors  were  associated  with  albedo  retrievals,  where  shaded 
facets  produced  albedo  estimates  up  to  15%  lower  than  actual  values.  Errors  increased  with  tilt  angles  up  to  20° 
but  decreased  with  sea  and  sky  turbidity.  Averaging  sunlit  and  shaded  pixels  before  running  the  inverse  model 
reduced  the  uncertainty  of  bathymetric  and  albedo  estimates  to  about  2  and  5%,  respectively,  comparable  to  previous 
field  evaluations  of  the  inversion  model. 


During  the  latter  half  of  May  and  the  beginning  of  June 
in  1998,  1999,  and  2000,  the  Coastal  Benthic  Optical  Prop¬ 
erties  (CoBOP)  study  occurred  in  the  waters  adjacent  to  Lee 
Stocking  Island  (LSI),  Bahamas.  One  of  the  goals  of  this 
study  was  to  provide  improved  measurements  and  under¬ 
standing  of  the  effect  of  different  bottom  types  on  remotely 
sensed  observations.  Remote  (aircraft,  satellite,  shipboard) 
and  in  situ  observations  of  the  water  and  bottom  were  made 
over  spatial  scales  that  ranged  from  a  few  centimeters  to 
kilometers.  As  part  of  the  interpretation  of  hyperspectral  im¬ 
agery  from  the  airborne  Portable  Hyperspectral  Imager  for 
Low-Light  Spectroscopy  (PHILLS)  sensor  (Davis  et  al. 
2002),  we  have  examined  the  contribution  of  sloping  or 
wave-scalloped  bottoms  (Carder  et  al.  2001)  on  upwelling 
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radiance  under  varying  conditions  of  illumination  and  tur¬ 
bidity. 

Recent  works  by  Lee  et  al.  (1999,  2001)  demonstrate  the 
utility  of  hyperspectral  radiance  data  collected  from  above 
the  sea  surface  for  determining  the  inherent  optical  proper¬ 
ties  and  chlorophyll  concentrations  of  the  water  column, 
even  for  shallow  waters.  This  inversion  method  (Lee  et  al. 
1999)  provides  estimates  of  bottom  depth  and  albedo  if  the 
assumption  is  made  that  the  bottom  is  horizontal  with  a 
Lambertian  bidirectional  reflectance  distribution  function 
(BRDF),  so  that  reflected  radiance  is  independent  of  the  in¬ 
cident  and  observational  angles.  For  shallow  waters  in  Tam¬ 
pa  Bay  (<5  m  depth),  bathymetric  estimates  accurate  to 
within  about  10%  of  measureiJ  depths  were  retrieved,  except 
for  regions  where  the  bottom  sloped  significantly  or  the  bot¬ 
tom-reflected  radiance  was  less  than  15%  of  the  total  mea¬ 
sured  radiance  (Lee  et  al.  2001). 

Zaneveld  and  Boss  (2003)  present  a  two-dimensional  an¬ 
alytical  model  of  upwelling  radiance  over  a  corrugated  bot¬ 
tom  to  examine  the  combination  of  bottom  morphology  and 
reflectance.  Although  their  model  does  not  include  skylight 
or  absorption  and  scattering  within  the  water  column,  it  does 
explain  some  of  the  differences  between  bottom  reflectance 
and  the  far  field  reflectance  when  the  measurement  is  made 
several  meters  off  the  bottom. 

Mobley  et  al.  (2003)  show  that  for  level,  ooid  sand  bot¬ 
toms  with  illumination  and  viewing  geometries  consistent 
with  PHILLS  aircraft  observations  at  LSI,  the  BRDF  is  es¬ 
sentially  Lambertian,  varying  only  about  4%  from  a  constant 
value  for  depths  greater  than  3  m  and  water  clarities  as  ob- 
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Fig.  1 .  PHILLS  images  of  (A)  Adderly  Cut  northwest  of  Lee  Stocking  Island  and  (B)  a  region 
west  of  Norman’s  Pond  Key  at  about  0915  h  EDT  on  17  May  2000.  RGB  bands:  663,  551,  468 
nm.  Transect  1  crosses  a  large,  shallow  sand  wave  feature  along  the  solar  plane.  Transects  2  and  3 
cross  sand  waves  along  and  orthogonal  to  the  solar  plane,  respectively. 


served  near  LSI.  Thus,  the  Lambertian  BRDF  assumption  of 
the  Lee  method  is  appropriate  for  PHILLS  observations  of 
these  sandy  areas.  Mobley  and  Sundman  (2003)  also  dem¬ 
onstrate  that  a  one-dimensional  radiative  transfer  model  with 
a  simple  correction  factor  (ID  X  CF)  for  solar  photons  can 
predict  up  welling  radiance  values  with  better  than  5%  ac¬ 
curacy  for  fiat,  gently  sloping  bottoms  (<10°  incline).  Their 
model  can  be  used  for  comparisons  with  backwards  Monte 
Carlo  (BMC)  modeled  radiance  reflected  from  a  corrugated, 
sawtooth  patterned  bottom  as  a  model  validation  step. 

Because  the  depth  gradient  of  a  bottom  can  be  estimated 
by  applying  the  Lee  inversion  method  to  hyperspectral  air¬ 
craft  imagery,  a  Mobley  and  Sundman  (2003)  correction  fac¬ 
tor  combined  with  the  Lee  inversion  method  should  improve 
accuracy  of  subsequent  bottom  albedo  estimates.  The  Mob¬ 
ley  correction  is  not,  however,  applicable  when  the  bottom 
is  not  flat  and  gently  sloping  or  when  the  bottom  is  corru¬ 
gated  within  the  area  viewed  by  an  aircraft  or  spacecraft 
pixel.  One  goal  of  this  paper  is  to  ascertain  the  accuracy  of 
water  depth  and  bottom  albedo  estimates  from  the  Lee  et  al. 
(2001)  inversion  method  for  corrugated  bottoms  in  shallow 
water. 

Bottom  reflectance,  water  depth,  optical  absorption  and 
scattering  all  contribute  to  the  variability  observed  in  scenes 
of  remotely  sensed  radiance.  PHILLS  aircraft  imagery  col¬ 
lected  during  the  CoBOP  study  at  LSI  show  repeating  wave¬ 
like  and  irregular  bottom  scour  features  (e.g.,  Fig.  1).  One 
explanation  for  the  wavelike  patterns  might  be  the  focusing 
of  sunlight  on  the  bottom  by  surface  waves,  as  discussed  by 


Zaneveld  et  al.  (2001).  However,  most  of  the  wavelike  pat¬ 
terns  apparent  in  the  PHILLS  imagery  have  wavelengths  >5 
m,  which  is  inconsistent  with  surface  waves  generated  in  a 
shallow,  short-fetch  region.  For  that  reason,  many  of  the 
wavelike  patterns  of  Figs.  1  and  2  are  thought  to  be  caused 
by  sand  waves,  and  the  presence  of  these  bedform  structures 
was  confirmed  by  diver  observations  and  were  documented 
with  microtopographic  images  of  the  bottom  measured  by 
the  Real-time  Ocean-Bottom  Optical  Topographer  (ROBOT, 
e.g.,  Fig.  3). 

In  addition  to  the  PHILLS  and  ROBOT  data,  a  variety  of 
in  situ  optical  measurements  of  the  water  column  and  bottom 
were  collected  near  LSI.  Although  the  spatial  resolution  of 
PHILLS  was  about  1.25  m,  the  resolution  of  the  ROBOT 
sensor  was  several  centimeters.  In  situ  vertical  profiles  of 
optical  properties  could  consist  of  a  series  of  centimetex- 
scale  measurements;  however,  these  profiles  or  ship-based 
measurements  of  remote  sensing  reflectance  might  have  been 
collected  kilometers  apart.  A  three-dimensional  BMC  model 
has  been  developed  to  combine  these  observations  in  a  nu¬ 
merical  simulation  to  determine  the  effects  of  bottom  topog¬ 
raphy  on  depth  and  albedo  estimates  using  hyperspectral  air¬ 
craft  data.  This  model  is  used  to  determine  whether  the  high 
contrast  of  the  wavelike  patterns  is  due  to  bottom  albedo 
differences  (e.g.,  heavy  minerals,  algae,  or  detritus  more 
concentrated  in  the  troughs  than  on  the  crests)  or  unequal 
illumination  (e.g.,  shadow  effects)  of  sand  waves. 

The  radiances  computed  by  the  BMC  model  over  a  uni¬ 
form  horizontal  bottom  are  comparable  to  radiances  pre- 
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Fig.  2.  PHILLS  remote  sensing  reflectance  transects  at  channels 
402,  559,  and  622  nm  over  sand  wave  and  algal  features  from  Fig. 
1  (Transect  1  from  Fig.  IB;  Transects  2,  3  from  Fig.  lA). 


dieted  using  Mobley’s  (1994)  Hydrolight  model.  The  BMC 
model  was  also  initiated  with  a  predominantly  oblique  light 
field,  fixed  average  water  depth,  and  a  corrugated  bottom  of 
uniform  albedo  to  mimic  the  environment  viewed  by  the 
PHILLS  sensor.  The  resulting  modeled  upwelling  radiances 
from  the  illuminated,  shaded,  or  composite  bottom  slopes 
were  then  used  as  surrogate  spectral  radiances  for  the  Lee 
reflectance  inversion  method.  The  resulting  estimates  of 
depth  and  albedo  were  compared  to  the  original  values  to 
evaluate  the  accuracy  of  the  Lee  inversion  method  when 
used  over  shallow,  corrugated  bottom  types.  The  BMC  mod¬ 
el  was  also  executed  using  a  sinusoidal  bottom  shape  for 
comparison  to  the  corrugated  and  flat  bottoms. 

Data 

Various  bottom  types  are  present  in  the  LSI  region,  in¬ 
cluding  ooid  sand,  coral  sand,  seagrass,  algal  mats,  and  ben¬ 
thic  diatoms.  Most  sand  albedos  or  bottom  reflectance  values 
at  LSI  ranged  between  0.2  and  0.6,  whereas  for  areas  with 
concentrations  of  either  detritus  or  algae,  the  reported  albe¬ 
dos  were  lower. 

PHILLS  is  a  hyperspectral  sensor  that  was  flown  over  LSI 
during  the  CoBOP  field  program  at  about  3,000  m  altitude, 
providing  a  spatial  pixel  resolution  of  about  1.25  m  (Davis 
et  al.  2002).  PHILLS  was  flown  when  the  solar  zenith  angles 
were  in  the  range  of  40-60°  in  order  to  avoid  sun  glint 


Fig.  3.  Three-dimensional  bathymetry  collected  from  the  RO¬ 
BOT  on  21  May  1999  over  an  ooid  sand  shoal  near  Adderly  Cut. 
Large  waves  have  a  sawtooth  pattern,  whereas  smaller  waves  are 
more  regular,  as  on  the  right. 


reflected  from  the  sea  surface.  Hence,  according  to  Snell’s 
Law,  the  subsurface  solar  illumination  was  directed  from 
about  32-40°  from  the  vertical.  Spectral  images  collected  by 
PHILLS  over  the  LSI  waters  include  areas  with  contrasting 
spatial  structures,  which  appear  to  be  caused  by  oblique 
lighting  of  sand  waves  on  the  bottom  (Fig.  1).  Contrast  is 
the  ratio  between  the  difference  and  the  average  of  radiance 
from  two  adjacent  areas.  In  addition  to  the  contrast  caused 
by  differing  illumination  of  sloping  wave  facets,  detritus  and 
algae  with  their  smaller  albedos  were  sometimes  found  on 
the  lee  sides  or  troughs  of  sand  waves,  modifying  the  con¬ 
trast  in  a  different  way.  Thus,  images  from  airborne  sensors 
such  as  PHILLS  provide  a  measure  of  the  effects  on  water¬ 
leaving  radiance  because  of  variability  in  both  the  topogra¬ 
phy  and  albedos  of  sand  bottoms  for  the  LSI  study  site. 

Transects  extracted  from  an  example  PHILLS  image  (Fig. 
1)  reveal  some  scour  patterns  (Fig.  2,  Transect  1)  and  wave¬ 
like  features  (Fig.  2,  Transects  2,  3).  Bottom  wave  facets, 
shaded  or  tilting  away  from  direct  sunlight,  appear  not  only 
darker  but  more  bluish  in  color  than  the  more  directly  illu¬ 
minated  facets.  For  example,  Transect  1  is  in  a  very  shallow 
region  and  the  blue  light  contrast  between  facets  was  smaller 
than  that  for  green  or  red  light,  probably  because  of  the 
reduction  of  direct  yellowish  solar  illumination  for  facets 
tilted  away  from  the  sun.  These  shaded  facets,  however,  are 
still  fully  illuminated  by  diffuse  blue  skylight.  Transect  2 
exhibits  a  similar  reduction  in  the  image  contrast  at  blue 
wavelengths. 

An  increase  in  the  contrast  exhibited  in  imagery  of  sand 
waves  can  also  be  caused  by  the  propensity  for  algae  to  grow 
on  the  more  quiescent,  lee  sides  of  sand  waves  in  high- 
current  areas  or  by  dark,  heavy  minerals  located  in  the 
troughs  of  sand  waves  (Fig.  1;  Fig.  2,  Transect  3).  In  these 
instances,  the  confrast  is  similar  for  blue  wavelengths  if  not 
greater  than  for  green  wavelengths.  This  would  be  expected 
because  algae  absorb  blue  light  more  completely  than  green 
light.  The  contrast  viewed  in  Transect  3  is  more  an  effect  of 
true  albedo  differences  than  of  interactions  between  illumi¬ 
nation  geometry  and  topographic  features  because  the  solar 
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plane  is  parallel  to  the  wave  crests,  equalizing  the  solar  ra¬ 
diance  incident  on  both  sides  of  the  crest  of  a  sand  wave. 

For  clean  waters  and  clear  skies,  scattered  light  near  the 
bottom  is  expected  to  be  predominantly  blue  because  of  the 
inverse-wavelength  (A""^)  behavior  of  Rayleigh  scattering  of* 
the  sky  and  water  molecules  and  the  reduced  absorption  of 
pure  water  at  blue  wavelengths  relative  to  longer  wave¬ 
lengths  (Pope  and  Fry  1997).  This  increases  the  apparent 
blue  reflectivity  of  indirectly  illuminated  bottoms  because 
this  scattered  light  is  similarly  reflected  from  both  sides  of 
a  wave,  whereas  yellow-rich,  direct  solar  light  is  reduced  on 
the  shaded  sides  of  sand  waves.  However,  heavy  minerals 
and  algae  absorb  as  much  or  more  light  at  440  nm  as  at  550 
nm.  The  blue-rich  reflectance  spectra  for  many  of  the  dark 
areas  of  wavelike  bottom  patterns  for  the  imagery  implies 
that  such  periodic  contrast  variations  were  not  caused  by 
increased  concentrations  of  benthic  diatoms  or  minerals  in 
the  troughs  of  the  waves  and  that  a  quantitative  evaluation 
is  needed  of  oblique  illumination  effects  on  a  rough  bottom. 

The  topographic  structure  representative  of  the  bottom  at 
LSI  is  required  to  parameterize  models  of  the  light  field.  The 
ROBOT  is  a  range-measuring  system  consisting  of  a  laser 
fan  beam  projected  onto  the  bottom  and  viewed  obliquely 
by  a  frame-processing  camera  as  described  in  Carder  et  al. 
(2001).  It  was  mounted  on  an  Ocean  Explorer— class  auton¬ 
omous  underwater  vehicle  provided  by  Florida  Atlantic  Uni¬ 
versity  (Smith  et  al.  1995)  and  nominally  operated  2  to  3  m 
above  the  bottom,  providing  a  bathymetric  swath  approxi¬ 
mately  1.5  to  2  m  wide. 

Some  bedforms  typical  of  the  ooid  sand  region  northwest 
of  LSI  (bright  area  containing  Transects  2  and  3,  Fig.  lA) 
are  presented  in  Fig.  3  as  measured  by  the  ROBOT.  The 
water  depth  in  this  area  is  about  4  m,  and  sand  waves  here 
appear  to  be  long-crested.  The  ROBOT  transect  shows  max¬ 
imum  bottom  slope  angles  of  <17°,  with  wavelengths  of  the 
sand  features  ranging  up  to  20  m.  Note  that  the  form  of  large 
bottom  features  might  not  be  sinusoidal,  but  sawtoothed  in 
shape,  whereas  the  small  wavelike  features  are  more  sinu¬ 
soidal. 

Discrete  water  samples  and  apparent  and  inherent  optical 
properties  were  measured  in  the  waters  surrounding  LSI.  The 
BMC  model  used  chlorophyll  concentration  as  an  input  pa¬ 
rameter,  with  all  optical  properties  covarying  with  it  as  in 
Morel  (1988).  Chlorophyll  a  (Chi  a)  concentrations  mea¬ 
sured  in  the  waters  immediately  adjacent  to  LSI  ranged  from 
O.IO  to  0.16  mg  m“^  with  an  average  value  about  0.14  mg 
m"\  We  evaluate  turbidity  effects  on  the  conti*ast  between 
brightly  and  dimly  illuminated  facets  of  sand  waves  due  to 
illumination  geometry  by  varying  Chi  a  concentrations  from 
0.01  to  0.5  mg  m-^  a  larger  range  than  observed  at  LSI,  in 
order  to  more  completely  evaluate  turbidity  effects. 

Models 

Mobley  et  al.  (2003)  discusses  the  necessity  of  using  a 
three-dimensional  radiance  transfer  (RT)  model  instead  of  a 
general  one-dimensional  (ID)  RT  model,  such  as  his  Hydro¬ 
light  model,  to  deal  with  the  problem  of  nonuniform  bot¬ 
toms.  They  also  discuss  a  simple  correction  factor  (ID  X 


CF)  that  increases  Hydrolight’s  predictive  accuracy  for  cases 
with  sloping  bottoms.  They  conclude  that,  for  a  bottom  slope 
angle  6^  <  10°,  a  ID  RT  model,  used  with  the  ID  X  CF 
correction  factor,  can  predict  the  nadir-viewed  radiance 
with  less  than  5%  uncertainty  for  solar  photons.  For  the  case 
of  a  corrugated  bottom  in  shallow  water,  however,  the  re¬ 
flected  radiance  between  adjacent  bottom  facets  plays  a  sig¬ 
nificant  role  in  as  do  internal  reflections  from  the  sea 
surface  and  diffuse  skylight.  Tliis  study  attempts  to  explore 
the  limitations  of  using  the  Mobley  ID  X  CF  approach  to 
predict  the  contrast  caused  by  the  sandy,  corrugated  bottoms 
when  sunlight  and  skylight  are  both  included.  We  compare 
the  Mobley  ID  X  CF  results  for  a  flat,  tilted  bottom  of  10° 
with  BMC  results  for  sawtooth-patterned  sand  waves  of  sim¬ 
ilar  slope  angles  to  explore  the  limitations  of  the  ID  X  CF 
approach. 

A  simplified  shallow-water  bottom  with  sand  waves  and 
the  associated  coordinates  are  illustrated  in  Fig.  4A.  The 
corrugated  bottom  is  modeled  as  continuous  wavelike  ramps 
repeating  in  the  x-direction  and  labeled  as  Ramp_2>  Ramp„„ 
Rampo,  Rampi,  and  Ramp2.  Each  simplified  wave  is  long- 
crested  and  orthogonal  to  the  solar  plane  (x-axis).  Because 
the  horizontal  variance  of  radiance  incident  on  the  sea  sur¬ 
face  is  negligible,  we  can  focus  on  a  single  ramp,  Rampo, 
and  assume  that  it  is  bounded  by  two  virtual  walls,  as  shown 
in  Fig.  4B.  If  a  photon  escapes  through  one  wall  at  coordi¬ 
nate  {x,  y,  z),  another  photon  traveling  in  the  same  direction 
can  be  supplied  immediately  to  the  same  position  (-x,  y,  z) 
on  the  opposite  wall  to  retain  the  photon  in  the  calculation 
cell.  The  optical  pathway  e  in  Fig.  4B  provides  an  example. 
This  assumption  does  not  change  the  computational  cost  of 
the  model.  It  is,  however,  convenient  for  the  program  code. 
Figure  4B  also  illustrates  the  optical  pathways  that  contribute 
to  the  sensor-detected  radiance  and  the  geometric  specifica¬ 
tions  of  the  bottom. 

The  BMC  model  is  an  extension  of  the  work  of  Liu 
(2000).  Using  radiance  reciprocity  arguments  (Gordon 
1985),  a  beam  of  light  is  simulated  by  a  very  large  number 
(bundle)  of  photons  emitted  from  the  sensor  to  the  water- 
body  (backwards).  A  set  of  bio-optical  models,  as  used  in 
Hydrolight,  is  employed  to  parameterize  the  inherent  optical 
properties  (lOPs)  as  a  function  of  the  Chi  a  concentration, 
with  the  water  column  assumed  to  be  homogeneous  (Case  1 
waters).  The  path  length  before  a  photon  hits  another  particle 
is  a  ftinction  of  the  beam  attenuation  coefficient  (Gordon 
1994).  Once  a  collision  occurs,  the  intensity  of  the  photon 
bundle  is  attenuated  by  multiplying  by  the  single-scattering 
albedo.  The  direction  of  bundle  travel  is  modified  by  scat¬ 
tering  and  the  ray  tracing  is  continued.  The  polar  scattering 
angles  relative  to  the  photon  direction  are  determined  by  the 
normalized  volume  scattering  phase  function  (Mobley  1994), 
and  the  azimuthal  scattering  angles  are  randomly  distributed 
(Gordon  1994).  Note  that  the  scattering  contribution  from 
pure  water  is  calculated  analytically,  whereas  the  contribu¬ 
tion  from  particulate  matter  is  specified  by  the  averaged  par¬ 
ticle  phase  function  (Mobley  1994).  If  the  photon  hits  the 
bottom  before  it  hits  another  particle,  the  position  of  inter¬ 
section,  as  well  as  the  bottom  slope,  is  recorded.  Tlie  bottom 
is  assumed  to  have  a  Lambertian  distribution  of  reflected 
radiance  with  a  bottom  albedo  or  irradiance  reflectance,  Pb. 


Illumination  effects  on  faceted  bottoms 


359 


Direct  solar  radiance 


Optical  pathways 
a:  hit  bottom  and  get  absorbed 
b:  hit  bottom  and  reflect 
c:  hit  particle  and  scatter 
d:  hit  particle  and  get  absorbed 
e:  hit  surface  and  reflect 
f:  hit  surface  and  refract 


Rampo 


Fig.  4.  (A)  Representation  of  a  shallow  bottom  with  sand  waves 

and  the  associated  coordinates,  with  d  defined  as  the  nadir  angle 
and  4>  as  the  azimuth  angle  and  with  the  x-axis  directed  toward  the 
sun.  (B)  Illustration  of  optical  pathways  that  contribute  to  the  sen¬ 
sor-detected  radiance  and  the  geometric  specification  of  bottom.  H 
is  the  mean  depth  of  the  water  column,  is  the  height  of  the 
ramp,  and  6^^  is  the  elevation  angle  of  the  ramp  from  the  hon- 
zontal.  Symmetry  allows  photons  leaving  the  right  wall  to  re-enter 
the  same  volume  on  the  left  (ray  e). 


Therefore,  the  photon  will  be  reflected  from  the  position  of 
intersection  into  any  direction  of  the  hemisphere  above  the 
local  bottom  facet,  and  its  intensity  is  attenuated  according 
to  the  irradiance  reflectance,  (Mobley  1994).  Ray  tracing 
continues,  with  special  attention  paid  to  whether  the  photon 
hits  an  adjacent  bottom  facet  before  it  hits  another  particle 
or  the  air-sea  interface.  The  sea  surface  is  assumed  flat.  If 
the  photon  hits  the  air-sea  interface,  Fresnel’s  formula  and 
the  laws  of  geometrical  optics  (Mobley  1 994)  are  used  to 
determine  the  photon’s  fate.  If  the  photon  is  refracted 
through  the  air-sea  surface,  the  remaining  intensity,  as  well 
as  its  direction,  is  registered.  If  the  photon  is  reflected  from 
the  air-sea  surface,  ray  tracing  is  continued. 

Above-surface  solar  irradiance  was  simulated  by  the  mod¬ 
el  of  Gregg  and  Carder  (1990),  and  the  sky  radiance  distri¬ 


bution  was  calculated  using  the  model  of  Harrison  and 
Coombes  (1988).  The  default  atmospheric  visibility  was  set 
to  15  km.  The  BMC  uses  the  same  sky  radiance  model  as 
Hydrolight  (Mobley  1994)  to  allow  direct  comparison  of  the 
BMC  and  Hydrolight  results.  The  sea  surface  is  assumed  flat, 
and  the  lOPs  of  the  water  column  are  assumed  homoge¬ 
neous.  A  set  of  bio-optical  models  for  Case  1  waters,  as  used 
in  Hydrolight,  is  employed  to  parameterize  the  lOPs  as  a 
frinction  of  the  Chi  a  concentration.  The  bottom  is  assumed 
to  have  a  Lambertian  distribution  of  reflected  radiance  with 
a  bottom  albedo  or  irradiance  reflectance,  pi,.  The  geometry 
of  the  ramp  is  specified  by  the  ramp  angle,  6^,,^,  and  the 
ramp  height,  The  radiance  field  can  then  be  computed 
at  any  location  in  the  water  column  for  any  viewing  angle. 

The  hyperspectral  inversion  method  (Lee  et  al.  1999, 
2001)  uses  a  priori  assumptions  about  water  column  con¬ 
stituents.  It  iteratively  varies  water  column  constituents, 
depth,  and  bottom  albedo  until  a  best  fit  is  obtained  to  the 
observed  hyperspectral  radiance  shape.  The  model  deter¬ 
mines  the  length  that  a  water  column  must  be  for  a  photon 
path  (sea  surface  to  bottom  to  surface)  to  filter  the  incident 
light  enough  to  provide  the  color  observed.  Because  of  the 
strong  absorption  of  water  in  the  red  wavelengths,  if  the 
color  of  water-leaving  radiance  is  red-rich,  the  water  column 
should  be  shallower  than  if  the  color  is  blue-rich. 

Because  the  Lee  inversion  technique  uses  an  optimization 
approach  to  best  fit  the  measured  remote  sensing  reflectance 
curves  with  model  curves,  minimizing  error  in  one  part  of 
the  spectrum  can  at  times  increase  it  elsewhere.  The  errors, 
then,  do  not  necessarily  scale  monotonically  with  wave  slope 
or  sun  angle.  For  the  present  study,  the  BMC  model  was 
constrained  to  six  spectral  channels.  Usually  the  inverse 
model  uses  at  least  18  spectral  channels  from  400  to  800 
nm. 


Results  and  discussion 

BMC  models  are  well  suited  to  the  task  of  evaluating  two- 
or  three-dimensional  environmental  simulations  (Gordon 
1985;  Reinersman  and  Carder  1995;  Reinersman  et  al.  1998; 
Mobley  et  al.  2003).  The  shape  of  recurring,  long-crested 
sand  waves  can  be  approximated  with  simple  ramp  pairs,  as 
shown  in  Fig.  4B. 

For  illustration.  Fig.  5  shows  the  BMC  model  results  for 
an  optically  clean,  10-m  water  column  where  the  corrugated 
bottom  has  a  ramp  height  of  0.8  m  and  ±10°  slopes.  The 
clean  water  is  an  extreme  case  because  diffuse  scattering  is 
minimized  along  the  path  of  a  photon.  Even  so,  the  diffuse 
contributions  from  the  sky  and  water  paths  are  significant, 
as  evidenced  by  the  reduced  400-nm  contrast  between  di¬ 
rectly  and  indirectly  illuminated  wave  facets.  Note  that  re¬ 
sults  for  ±10°,  tilted  flat  bottoms  (Hydrolight  with  ID  X 
CF;  Mobley  and  Sundman  2003)  where  the  total  incoming 
irradiance  to  the  sea  surface  is  the  same  as  for  the  BMC 
model  case,  nearly  match  the  corrugated  BMC  results  for 
individual  facets  at  560  nm.  At  400  nm,  however,  the  tilted, 
flat  bottom  provides  greater  contrast  than  the  BMC  corru¬ 
gated  results.  The  BMC  predictions  for  a  10°  corrugated  bot¬ 
tom  also  produces  significantly  less  contrast  at  400  nm  than 
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Fig.  5.  (A)  Green  and  (B)  blue  radiance  modeled  for  the  sub¬ 

surface  view  above  a  sand  wave  crest  (Fig.  4)  with  //  -  10  ^rniup 
=  0.8  m,  Pb  =  50%,  and  extremely  clear  water  {see  case  I  in  Morel 
[1988]:  Chi  a  =  0.01  mg  m"^).  The  sand  wave  crests  are  oriented 
perpendicular  to  the  solar  plane,  and  the  solar  zenith  angle  is  60°. 
Solid  squares  are  model  results  with  conugated  ramp  bottom;  open 
circles  are  results  for  a  sinusoidal  bottom.  Dashed  lines  are  Hydro¬ 
light  results  with  correction  factors  for  flat  bottoms  (solid  line) 
sloped  at  10°  toward  (upper  line)  and  away  (lower  line)  from  the 
sun. 


at  560  nm  because  of  shadow  infilling  by  diffuse  (mostly 
Rayleigh)  light.  Recall  that  the  ID  X  CF  modification  to 
Hydrolight  applies  only  to  direct  photons  and  does  not  in¬ 
clude  skylight,  which  represents  57%  of  the  incoming  light 
at  400  nm  and  41%  of  the  light  at  560  nm.  For  a  visibility 
factor  of  50  km,  these  ratios  become  42  and  20%,  respec¬ 
tively,  which  would  increase  the  contrast  of  the  sand  wave 
facets.  The  up  welling  radiance  can  be  shown  to  be  maximal 
on  the  illuminated  side  of  the  corrugation  crest,  where  the 
incidence  angle  equals  the  reflectance  angle  for  rays  hitting 
the  bottom.  With  a  solar  zenith  angle  of  60°,  the  subsurface 
solar  zenith  angle  is  about  40°,  so  the  bottom  slope  project¬ 
ing  the  brightest  return  toward  zenith  is  about  20°  from  hor¬ 
izontal. 

Note  in  Fig.  5  that  the  10°  sawtooth  pattern  produces  a 
smaller  maximal  reflected  signal  than  does  the  sinusoidal 
pattern,  which  has  a  maximum  bottom  slope  angle  of  14.1°^ 
which  is  closer  to  the  brightest  bottom  slope  angle  of  20° 


Table  1 .  List  of  the  computational  conditions  used  in  the  back¬ 
wards  Monte  Carlo  model  for  the  puiposes  of  understanding  how 
the  contrast  is  influenced  by  various  factors.  Each  factor  of  the 
standard  case  is  varied  one  at  a  time  as  the  comparative  case. 


Standard 

case 

Comparative  case 

Solar  zenith  angle,  (degrees) 

60 

0,  30,  60 

Depth  of  water  column,  H{m) 

5 

2,  5,  10 

Chi  a  concentration  (mg  m"^) 

0.0 1 

0.01,  0.1,  0.5 

Lambertian  bottom  reflectance, 

0,5 

0.2,  0.5 

Atmospheric  visibility,  (km) 

15 

5,  15,  50 

Ramp  angle,  (degrees) 

10 

10,  15,  20 

than  is  10°.  BMC  model  simulations  for  bottom  slope  angles 
up  to  20°  are  made  in  order  to  include  angles  providing 
maximal  radiance  from  the  bottom  to  a  nadir-viewing  aircraft 
sensor.  From  Fig.  5,  it  is  clear  that  the  effect  of  the  water 
depth  variation  of  0.8  m  (slope  from  viewing  angle  0-24°) 
has  a  smaller  effect  on  the  upwelling  radiance  than  has  the 
bottom  slope  angle  (slope  from  viewing  angle  +0°  to  -0°), 
especially  at  400  nm.  Note  for  the  sawtooth  sand  waves,  the 
crest  versus  trough  values  for  the  sunlit  side  (viewing  angle 
0-24°,  Fig,  5)  decrease  because  of  a  range  increase  from  sun 
to  bottom  to  sensor  of  roughly  1.9  m,  including  slant-path 
effects.  The  additional  absorption  at  560  nm  along  that  path 
difference  is  about  12%,  which  is  very  close  to  the  13% 
decrease  observed  from  0  to  24°  in  Fig.  5.  The  shade  effect, 
on  the  other  hand,  is  easy  to  observe  by  comparing  LJ0°) 
to  LJQV)  because  both  are  at  the  same  depth.  This  decrease 
is  about  17%  at  560  nm,  or  nearly  50%  larger  than  the  depth 
effect  for  the  sunlit  facet. 

The  effect  of  depth  difference  due  to  sand  waves  has  only 
a  3%  effect  at  400  nm  because  of  tlie  small  absorption  co¬ 
efficient  (—1%)  for  these  clear  waters.  The  shade  effect  is 
11%,  or  about  two  thirds  of  the  effect  seen  at  560  nm  be¬ 
cause  of  the  shadow  infilling  by  additional  diffuse  sky  and 
water  path  radiance  at  blue  wavelengths. 

To  understand  how  the  contrast  between  adjacent  wave 
facets  is  influenced  by  various  factors,  a  series  of  simulations 
(six  spectral  radiances)  were  executed  based  on  the  com¬ 
putational  conditions  listed  in  Table  1 .  All  factors  of  a  stan¬ 
dard  case  are  varied  one  at  a  time  for  sand  waves  with  a 
sawtooth  pattern.  A  hypothetical  sensor  is  placed  above  a 
corrugation  crest,  just  below  the  air-water  surface  (depth  z 
=  0"),  providing  a  nadir  view  (0J  of  either  the  sunny  {6^  = 
0+)  or  shady  sides  {6^  =  0")  of  the  crest.  The  contrast  is 
defined  as  the  ratio  between  the  difference  and  the  average 
of  radiance  from  both  sides  near  the  crest,  as  detected  by 
the  hypothetical  sensor. 

,  l^(^v  =  0^)  -  =  0~)l  Hi 

contrast  (a  =  0  )  =  o.5[I(0„  =  O')  +  =  0^)] 

Spectral  contrast  changes  due  to  environmental  variabil¬ 
ity — Figure  6A  shows  that  the  solar  angle  significantly  af¬ 
fects  the  contrast  of  illuminated  versus  shaded  sides  of  10° 
corrugations.  Note  that  with  the  sun  at  zenith,  the  corruga- 
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Fig.  6.  Spectral  contrast  variations  for  sawtooth  sand  waves 
with  environmental  conditions  listed  in  Table  1 .  Effects  due  to  (A) 
solar  zenith  angle,  6s^„;  (B)  water  column  depth,  H;  (C)  Chi  a  con¬ 
centration  and  covarying  lOPs;  (D)  bottom  albedo,  (E)  atmo¬ 
spheric  visibility  at  550  nm,  V^;  and  (F)  ramp  angle  of  wave  facets, 

^romp* 


tions  would  be  difficult  to  discern  because  their  contrast 
would  be  near  zero. 

For  water  less  than  10  m  deep,  the  effect  of  bottom  depth 
on  corrugation  contrast  is  relatively  minor  (Fig.  6B).  For 
deeper  clear  waters,  however,  the  contrast  decreases  slightly, 
especially  at  the  Rayleigh-rich  blue  wavelengths.  This  is  be¬ 
cause  of  the  increase  in  diffuse  molecular  scattering  arising 
from  the  longer  light  paths  of  deeper  waters.  Similarly,  as 
water  turbidity  increases,  contrast  decreases.  In  this  BMC 
model,  increasing  turbidity  is  caused  by  an  increase  in  chlo¬ 
rophyll  concentration  (Fig.  6C). 

Increasing  the  bottom  albedo  results  in  a  slight  decrease 
in  contrast  (Fig.  6D).  This  is  because  both  the  shaded  and 
illuminated  areas  are  subject  to  multiple  scattering  from  the 
water  column,  to  reflected  radiance  from  an  adjacent  facet 
and  to  internal  reflectance  from  the  sea  surface. 

Decreasing  the  atmospheric  visibility  (transparency)  from 
50  to  5  km  (Fig.  6E)  makes  the  light  field  above  the  sea 
surface  more  diffuse.  The  BMC  model  predicts  a  reduction 
in  contrast,  as  one  would  expect. 

Increasing  the  facet  slope  of  bottom  corrugations  where 
they  are  obliquely  illuminated  significantly  increases  the 
contrast  (Fig.  6F).  With  the  sun  at  a  large  solar  zenith  angle 
(60°)  and  shining  through  clear  water  and  sky,  the  contrast 
between  facets  of  sand  waves  can  become  as  large  as  45% 


for  bottom  facet  slopes  of  20°.  Recall  that  a  nadir  view  of  a 
bottom  facet  with  a  20°  slope  toward  the  sun  is  optimal  for 
maximal  reflectance  for  subsurface  illumination  at  a  40°  so¬ 
lar  zenith  angle. 

Spectral  bathymetric  and  albedo  retrievals — The  radianc¬ 
es  generated  from  the  BMC  simulations  discussed  above 
were  subsequently  used  as  inputs  for  the  inversion  method 
of  Lee  et  al.  (1999,  2001).  This  inversion  of  the  Lee  et  al. 
(1998)  remote  sensing  reflectance  modef  provides  estimates 
of  water  lOPs,  water  depth,  and  bottom  albedo.  Subsequent¬ 
ly,  these  estimated  depths  and  albedos  are  compared  to  the 
original  values  used  by  the  BMC  model.  The  estimates  are 
not  expected  to  exactly  match  the  original  values  because 
the  Lee  inversion  assumptions  include  a  horizontal  bottom, 
whereas  the  BMC  model  has  produced  near-nadir  views  of 
the  illuminated  and  shaded  bottom  corrugations. 

Table  2  presents  errors  in  bathymetric  and  albedo  esti¬ 
mates  for  the  BMC  simulations.  Error  percentages  are  pre¬ 
sented  as  variations  in  the  height  of  the  water  column,  H, 
and  bottom  albedo,  p,  for  the  sunny  and  shady  sides  of  the 
corrugation  crests.  These  simulations  were  for  long  sand 
waves,  (wavelength  >  10  m),  as  illustrated  by  the  ROBOT 
data  in  Fig.  3.  Errors  for  smaller  sand  waves  with  less  ver¬ 
tical  relief  and  slopes  will  be  smaller.  Also  depicted  in  Table 
2,  Case  1 ,  are  the  errors  obtained  if,  before  retrieving  model 
inversion  values,  the  illuminated  and  shaded  BMC  facets  are 
averaged  to  simulate  radiance  data  from  a  sensor  with  lower 
spatial  resolution. 

With  clear  water  and  a  corrugated  bottom,  the  Lee  model 
underestimates  water  column  depths  for  the  sunlit  facets  and 
overestimates  depths  for  shady  ones.  Most  errors  are  less 
than  5%,  except  where  slopes  approach  20°.  Averaging  the 
contrasting  radiances  provides  depths  accurate  to  within  3% 
for  all  cases  in  Table  2,  Case  1,  even  for  slopes  of  20°. 

The  error  in  the  bottom  albedo  can  reach  15%  for  shaded 
facets.  However,  except  for  very  shallow  water  depths,  ac¬ 
curacies  improve  to  about  5%  by  averaging  adjacent  con¬ 
trasting  radiances. 

Increasing  water  turbidity,  depth,  and  bottom  absoiption 
all  increase  the  relative  contribution  of  the  water-reflected 
photons  relative  to  bottom-reflected  photons,  decreasing  the 
influence  of  the  bottom  reflection  on  the  above-water  spec¬ 
tra.  Increased  water  turbidity  reduced  the  accuracy  of  albedo 
estimates  for  shaded  facets;  however,  it  increased  bathymet¬ 
ric  accuracies  for  both  sides  of  a  sand  wave  because  diffiise 
light  reduces  the  color  distortions  of  the  two  facets  from  that 
of  a  horizontal  bottom.  Turbidity  effects  should  not  be  dom¬ 
inant  error  sources  for  the  PHILLS  imagery  evaluated  here. 

Table  2  also  provides  results  for  a  hypothetical  water  type, 
where  the  light  absorption  by  dissolved  material  is  greater 
than  that  of  chlorophyll  at  440  nm  (Case  2  waters).  This  is 
more  representative  of  the  water  conditions  for  the  Adderly 
Cut  region  near  LSI  than  the  conditions  used  in  Table  2, 
Case  1 .  According  to  these  results,  bathymetry  errors  ca. 
<3%  and  albedo  errors  <15%  can  be  attributed  to  10°  sand 
waves  if  the  method  of  Lee  et  al,  (2001)  is  applied  to  ac¬ 
curate  spectral  images  of  the  Adderly  Cut  region.  For  spa¬ 
tially  averaged  data,  including  both  sunny  and  shaded  sides 
of  the  waves,  the  sand  wave  errors  reduce  to  less  than  1.4 
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Table  2  Error  percentages  of  bathymetric  and  albedo  retrievals  for  a  wide  variety  of  lighting  and  environmental  conditions  for  long 
waves,  such  as  illustrated  by  ROBOT  data  in  Fig.  3.  Case  1  waters,  of  which  the  lOPs  are  a  fonction  of  the  Chi  a  concentration  [Chi  o] 
based  on  a  set  of  bio-optical  models,  are  as  used  in  Hydrolight.  Gelbstoff-rich  Case  2  waters  have  an  extra  CDOM  component  added.  It 
only  contributes  to  absorption  and  not  to  scattering.  The  absorption  by  CDOM  is  set  to  be  2.5  times  that  of  the  chlorophyll  absorption  at 
a  reference  wavelength  of  440  nm. 


Error  of  bathymetric  and  albedo  retrievals  {%) 

^Sun  ^  [Chi  g] _ ^fainp _ Pb _ ^shacly _ _ ^average _ Pshady _ Psunny 


Case  1 

0 

30 

5 

0.0 1 

10 

60 

60 

2 

5 

0.0 1 

10 

10 

60 

0.01 

5 

0.1 

10 

0.5 

60 

10 

5 

0.01 

15 

20 

60 

5 

0.01 

10 

Case  2 

60 

2 

5 

0.1 

10 

10 

60 

2 

5 

0.1 

10 

10 

0.94 

1,28 

1.07 

0.5 

3.56 

‘  -1.58 

0.77 

5.86 

-2.69 

-0.16 

6.05 

-3.22 

-0.94 

0.5 

5.86 

-2,69 

-0.16 

5.74 

-3.68 

0.46 

5.86 

-2.69 

-0.16 

0.5 

6.26 

-1.44 

2.89 

1.35 

-0.39 

2.30 

5.86 

-2.69 

-0.16 

0.5 

9.02 

-3.82 

0.93 

11.78 

-4.89 

1.33 

0.2 

1.64 

-4.18 

-0.57 

0.5 

5.86 

-2.69 

-0.16 

-0.14 

-2.50 

-1.38 

0.2 

1.03 

-0.52 

0.28 

-3.03 

0.80 

-0.94 

2.20 

-2.81 

-0.42 

0.5 

-0.23 

-0.91 

0.48 

-2.05 

-0.69 

-1.17 

4.07 

7.77 

1.38 

6.76 

-5.77 

1.97 

-9.79 

2.07 

-5.77 

1.97 

-2.32 

3.29 

-5.77 

1.97 

-6.03 

2.85 

-12,00 

0.53 

-5.77 

1.97 

-10.03 

4.54 

-15.60 

3.87 

-13.49 

0.54 

-5.77 

1.97 

-15.29 

2.25 

-12.22 

3.98 

-14.89 

6.58 

-15.75 

-2.88 

-14.81 

2.83 

-14.43 

2.68 

5.88 

4.00 

-2.39 

-3.50 

-2.39 

I.IO 

-2.39 

1.38 

-1.87 

-2.39 

-2.04 

-4.28 

-5.91 

-2.39 


-6.42 

-3.98 

-4.27 

-9.12 

-3.62 

-5.55 


and  9.1%,  respectively.  These  errors  would  be  in  addition  to 
those  implicit  in  the  inversion  method  (Lee  et  al.  1999). 
Increasing  the  wave  slopes  will  increase  the  errors,  but  re¬ 
ducing  the  solar  zenith  angle  will  decrease  the  errors. 

Lee  et  al.  (1999)  found  that  their  estimates  of  depths  av¬ 
eraged  within  10%  of  the  measured  depths  when  using  data 
from  the  Florida  and  the  Bahamas  shelves.  For  Tampa  Bay 
in  Florida  (Lee  et  al.  2001),  bathymetric  accuracies  averaged 
about  10-12%,  except  where  the  bottom  slope  was  signifi¬ 
cant  or  the  bottom  signal  contributed  less  than  15%  to  the 
water-leaving  radiance.  A  slope  of  20®  inducing  additional 
errors  up  to  1 2%  in  bathymetry  and  1 6%  in  albedo  for  the 
shady  sides  of  waves  or  sloping  bottoms,  as  found  in  this 
study,  is  consistent  with  the  field  application  of  Lee  et  al. 
(2001). 

PHILLS  aircraft  imagery  of  the  shallow  waters  near  LSI 
contain  long-crested,  wavelike  bedforms  that  might  have  re¬ 
sulted  firom  shading  because  of  oblique  illumination  of  sand 
waves  or  decreased  albedo  in  the  trough  regions  of  the 
waves.  The  profiles  of  the  reflectance  patterns  observed  in 
the  PHILLS  imagery  (Fig.  2  transects)  are  angular  rather 
than  smooth,  suggesting  the  bottom  waveform  is  sawtoothed 
rather  than  sinusoidal.  The  BMC  model  uses  the  assumption 
that  the  bottom  is  formed  of  repetitive  symmetric  ramps, 
whereas  the  actual  bottom  profiles  measured  with  the  RO¬ 
BOT  appear  to  be  composed  of  as3mimetric  ramps. 

With  oblique  illumination,  the  model  provides  shade  ef¬ 
fects  consistent  with  patterns  observed  in  PHILLS  imagery 
along  Transects  1  and  2.  Notably,  the  contrast  at  blue  wave¬ 


lengths  was  lower  than  at  longer  wavelengths,  and  the  con¬ 
trast  ranges  were  generally  within  the  ranges  observed  in 
model  results.  When  the  BMC  model  was  run  with  a  60® 
solar  zenith  angle,  then  ramp  slopes  as  high  as  20®  and  at¬ 
mospheric  visibilities  as  high  as  50  km  were  required  to 
approximate  the  reflectance  conti'ast  observed  at  blue  wave¬ 
lengths  along  Transect  1.  However,  PHILLS  flew  when  the 
solar  zenith  angle  was  52®,  so  the  BMC  overestimate  of  con¬ 
trast  for  this  PHILLS  image  is  plausible.  The  observed  con¬ 
trasts  of  red  and  green  reflectance  were  much  larger  than 
predicted  by  the  BMC  model,  however,  so  the  height  of  the 
sand  waves  might  be  greater  than  the  model-assumed  height 
or  the  albedo  of  the  troughs  might  be  darker  than  those  of 
the  crests. 

The  profiles  of  Transect  3,  however,  are  inconsistent  with 
simple  topographic  forcing  of  the  obliquely  lit  light  field 
because  the  transect  is  normal  to  the  solar  plane,  minimizing 
topographic  effects,  and  because  blue  contrast  features  equal 
or  exceed  variations  observed  at  green  and  red  wavelengths. 
Variations  in  bottom  albedo  from  algal  accumulations  in 
sand  wave  troughs  would  provide  a  logical  explanation  for 
these  contrast  variations. 

Bottom  slope  effects  on  bathymetric  and  albedo  retrievals 
can  be  minimized  by  making  flight  operations  at  smaller 
solar  zenith  angles  (e.g.,  40®)  under  conditions  where  sand 
wave  crests  run  parallel  with  the  solar  principal  plane.  Fur¬ 
thermore,  if  the  slope  of  the  bottom  can  be  infen'ed  from 
sequential  bathymetric  retrievals  across  wave  forms  in  an 
image,  there  is  a  fair  chance  that  the  slope  illumination  ef- 
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fects  on  the  albedo  of  the  bottom  can  be  corrected,  at  least 
in  part  (e.g.,  Mobley  and  Sundman  2002).  However,  such  a 
research  effort  is  beyond  the  scope  of  this  paper  and  is  left 
for  future  consideration. 
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